Functional integrals for QCD at nonzero chemical potential and zero density 
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In a Euclidean space functional integral treatment of the free energy of QCD, a chemical potential 
enters only through the functional determinant of the Dirac operator which for any flavor is Ip 
+ m — /i/7o (where [if is the chemical potential for the given flavor). Any nonzero /i alters all of the 
eigenvalues of the Dirac operator relative to the fi = value, leading to a naive expectation that 
the determinant is altered and which thereby alters the free energy. Phenomenologically, this does 
not occur at T = for sufficiently small /i, in contradiction to this naive expectation. The problem 
of how to understand this phenomenological behavior in terms functional integrals is solved for the 
case of an isospin chemical through the study of the spectrum of the operator 70 (P+ m). The case 
of the baryon chemical potential is briefly discussed. 
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Consider QCD with two degenerate light flavors at zero temperature and a finite isospin chemical potential. At a 
phenomenological level much is known about this system The behavior at low chemical potential is of interest 
here. The energy density and isospin density are known to be zero for < m ff . As is increased above a 
second order phase transition to a pion condensed state is reached in which the energy density and isospin density 
both continuously increase from zero. It is pretty clear how to interpret these results in terms of eigenstates of the 
QCD hamiltonian. The critical value of hi corresponds to the state in the system with the smallest energy per unit 
isospin (namely a pion at rest). Unfortunately, there is no simple way to obtain these eigenstates starting directly 
from the QCD. Here the focus will be on understanding the phenomenological behavior in terms of a Euclidean space 
functional integral formulation of the theory. 

For simplicity of presentation, the discussion will be restricted to two-flavored QCD with degenerate flavors. The 
inclusion of additional heavy flavors is straightforward and does not alter any of the results. As will be clear, it is 
necessary to work at finite temperature and then study the limit as T — > 00 at a later stage. The free energy is 
given by Gi(T, = E — TS — ^7/3 where E, T, S and /ij are the energy, temperature, entropy density, isospin 
and chemical potential respectively. One can envision working in a finite but large box of volume V with appropriate 
boundary conditions. The thermodynamic limit V — > 00 can be taken at the end of the calculation and results may 
then be expressed in terms of intensive variables such £ = E/V, Q = G/V or pj = I3/V. The grand partition 
function, Zj(T/ij) = e ~/ 3G ( T .M/) ; w ith = 1/T can be represented as a functional integral, 

(1) 

where Sym is the Yang-Mills action and the functional integral is evaluated over a four dimensional box with a 
temporal length of (3 and with periodic boundary conditions for the gluons. The functional determinant is computed 
over the same box but with anti-periodic boundary conditions for the fermions and is taken to be over a single flavor. 
The absolute value squared arises since the functional determinant of the up quarks is the complex conjugate of the 
functional determinant of the down quarks provided the chemical potentials for the two flavors are equal and opposite 
H 0. Thus the only way that the presence of the chemical potential is felt is via the functional determinant of 
the Dirac operator. The determinant may be represented as the product of the eigenvalues of the Dirac operator 
det(#>+m-^ ) =H j X j where (#>+ m - ^ ) ipj = Xjipj . 

At this stage one is confronted by something of a puzzle. A priori, one expects that the eigenspectrum of the Dirac 
operator with fij — (and some fixed gauge configuration) is completely different from the spectrum with any nonzero 
fij in the sense that every eigenvalue is expected to be different; there is no mathematical reason for the eigenvalues 
not to depend on \i. In the absence of some conspiracy between the eigenvalues, one would therefore expect that the 
functional determinant with any nonzero fii would differ from the functional determinant [ii — and that this would 
occur for every gauge configuration. In turn, eq. JQ) , would lead one to expect that if the functional determinants all 
differ from their tt/ = counterparts, that Zi(T,/j,i) also would vary. Moreover, this expectation would seem to hold 
at any temperature including T = 0. Clearly, this expectation is wrong. At T — we know phenomenologically that 
the free energy is precisely its vacuum value for nonzero values of jii provided that < m^. The puzzle is simply 



Zi(T tl a) = / d[A] det(>- 
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how can we understand the vanishing physical effects of < in the context of a functional integral treatment. 

In one sense, this problem might be dismissed as, quite literally, much ado about nothing; the problem amounts 
to getting a mathematical understanding of nothing happening at the physical level. However, to understand the 
phase transition where something does occur, it is necessary to understand just how it is that nothing occurs below 
the transition. Accordingly, it seems fitting to name this the isospin "Silver Blaze" problem after the Arthur Conan 
Doyle story of that name in which the "curious incident" of a dog doing nothing in the night time provides Holmes 
with an essential clue. There are analogous Silver Blaze problems associated with the baryon chemical potential, the 
strangeness chemical potential or various linear combinations of these chemical potentials. 

The key to the solution of the isospin Silver Blaze problem is the study of eigenstates of 70 times the usual Dirac 
operator. The eigenspectrum of the usual Dirac operator, p+ m, has received considerable attention in both random 
matrix models)^ and lattice studies[{|; the intense interest is because this eigenspectrum provides essential insight into 
chiral symmetry breaking[6j. However 70 times the Dirac operator has not received such interest. There have been 
some studies of the so-called propagator matrix, which is a lattice version of 70 times the Dirac operator in connection 
with the Glasgow method^] and in a random matrix form||. Such treatments are not completely illuminating for the 
present purposes since the key result is manifest only in the continuum limit (or strictly speaking where the continuum 
limit has been taken at least in the time direction). In this limit, the eigenspectrum of 70 times the Dirac operator 
provides insight into the physics of pion condensation in much the same way as the eigenspectrum of the usual Dirac 
operator gives insight into chiral symmetry breaking. To see how this comes about consider the partition function 
of QCD in the presence of an isospin chemical potential as given in eq. with /zj solely affecting the functional 
determinant. It is a simple exercise in linear algebra to show that functional determinant can be represented as 



det($>+TO-^J =det(p+m) e LV / J (2) 

which shows that all of the effects of the chemical potential can be represented as the exponential of an integral of the 
trace of 70 times the usual Dirac operator. This trace can be evaluated via a sum over eigenvalues with the functions 
subject to anti-periodic boundary conditions. 

There are several features of the operator (jo(p+ m) — which greatly facilitate the computation of the trace. 
One useful property concerns the paring of eigenvalues. It is easily demonstrated that if A is an eigenvalue of 
("fo(fl+ m)), then so is —A*. Next notice that the time derivative term in the operator is not multiplying any 7 
matrices since 7q = 1. Moreover, also does not multiply a 7 matrix but simply enters as an additive constant and 
consequently has no role in determining the eigenfunctions; it merely shifts the eigenvalues. Together these imply that 

if (j (P+ m)) = A|V>), then (jo(p+ m) - 4) e Mt \ip) = e iuJt \ip). This means that for any given 

eigenfunction, another can be generated with an eigenvalue shifted by iuj by simply multiplying by a time-dependent 
phase. Note at this stage it is critical that the continuum expression is used for the time derivatives. Moreover, these 
phase factor cannot be arbitrary. Recall that the eigenfunctions are constrained to be anti-periodic in time with a 
periodicity of [3. Thus, the eigenfunctions of (^jo(]p+ m) — can be grouped into families denoted by two indices 
an index j representing an "intrinsic" eigenstate and an index n representing the phase factor satisfying: 



, , u't . ( Oj (2n + l)ir \ 
|Vin+i) = eT~ |^n> A jn = e J-Y +l (j3 + jg ) (3) 



where A jn , the eigenvalue of the operator, is broken up into a real and an imaginary part. This decomposition is 
unique if we impose the condition that — tt < tfij < ir. The form of eq. 10 is highly suggestive. Indeed, if one ignores 
the fact that <$>j is generically nonzero, it is identical to the form for the case of free noninteracting fcrmions; in the 
noninteracting case, €j simply represents the energy of a mode. Here, e^, the real part of eigenstates of 70 times the 
usual Euclidean Dirac operator, may be considered as a quasi-energy. Of course, the quasi-energies differ from the 
energies in the free fermion case in critical ways; they depend on the gauge field configuration and are not known 
analytically. 

The trace can be evaluated by first summing over the n for the eigenvalues (which is analogous to a usual Matsubura 
sum0 but here done for fields in the presence of a time-dependent background); a sum over the intrinsic j with ej > 
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is done subsuently: 



tr 



7o {Jp+m) 







tanh 
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Recall that the problem of interest is at zero temperature which corresponds to (3 — > oo. As this limit is approached, 
the asymptotic behavior of the hyperbolic tangent can be used to simplify the result of eq. Q : 



tr 



7o(^+m)-^ 



= £ 9{e s ) sign(^) {-pe(\^\ - 2 £j ) + i2cf> j 6(\^\ - 2e j )] + 0(e^ A ) 
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where A represents a characteristic hadronic scale. Note that this becomes exponentially small as j3 — ► oo and will 
have a negligible effect. Inserting eq. into eq. yields 



dct (p+ m-Hqs.^ 
det (p+m) 



= exp I -iY,^6{e j )Q{\^ I \-2e j ) I exp I § £ flfo) 0(M - 2e,) - 2 Cj ) + 0(e"^ A ) 



(6) 

In using eq. © rather than the full expression in eq. Q and neglecting the exponential, one is implicitly assuming 
that the zero temperature limit for the physical quantity of interest, the free energy, has a smooth behavior as the 
zero temperature limit is approached. This is highly plausible at finite chemical poitential but has not been directly 
proved from QCD. 

Equation © provides an essential ingredient for the resolution of the isospin Silver Blaze problem. The theta 
functions ensure that for any given gauge configuration in the QCD functional integral, as T = is approached, 
the functional determinant is unchanged from its /i/ = value, unless is greater than the quasi-energy of the 
minimum positive quasi-energy mode. Thus, eq. JfjJ explains how at zero temperature, functional determinants can 
remain unchanged despite all of the eigenvalues of the Dirac operator changing which is the core of the Silver Blaze 
problem. 

Of course, the value of the minimum positive quasi-energy depends on the gauge configuration. A full resolution of 
the isospin Silver Blaze problem depends on showing that nothing physical occurs at zero temperature unless fii > m n . 
Clearly this is the case, provided that the minimum positive quasi-energy is greater than or equal to for those 
configurations which contribute to the functional integral at zero temperature. A formal statement of this condition 
can be given in terms of the spectral density (as a function of e ) averaged over the gauge conditions. We denote this 
spectral density as p(e) 



f>(e) = '£S(e 
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(70GP+ m) - e)- 1 + ((-#+ m) 7o - e)' 1 + O^) 



(7) 



where ej is the j th quasi-energy for a given configuration, and the second equality follows from eq. (JSJ). We denote 
averaging over gauge configurations with the notation: 



(6) 



1 



d[A]6 detfjp 
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where the functional integral is evaluated with standard boundary conditions. The minimum relevant positive quasi- 
energy, e m i n , can be defined in the following way: (/3(e))o,o = if and only if |e| < e m i n . The Silver Blaze problem is 
then fully resolved if two conditions are met: 

i) (p( e ))o,fu = (p( e ))o,o for all < 2e mi „ 



To verify that these conditions are sufficient, use the definition of the free energy in terms of Zj and eqs. i ffl l) and Q 

to note that — — 2/ 2 de (/3(e))o,^j. Thus if T = 0, < and the two conditions are satisfied, then free 
energy is independent of /ii and the expectation value of the isospin vanishes . 
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It is straightforward to demonstrate that the two conditions are in fact met in QCD, provided that one assumes 
that no first order phase transition occurs for T = and \pi\ < m,. This assumption is innocuous — it is highly 
plausible a priori and is known not to occur in nature. To demonstrate the validity of i) use eqs. (J7J) and (JSJ to show 
that 



d_ 




— m)o^ = (p(c)[ I de'p(e')\ ) (9) 



Note that, the operator J Q 2 de' j o(e') is manifestly nonnegative as is the measure in eq. (JSJ). In general, if an operator 
A is manifestly nonnegative and it is averaged over a nonnegative measure and if (A) = it follows that (BA) = 0. 
Equation 10 and this general condition together imply that g^-(p(e))o,^j = if (L 2 de'p(e')) =0. This in turn 

implies the validity of condition i) provided that (p(e))o,/jj does not change discontinuously at some value of [ii (i.e., 
there is no first order transition). 

To establish conditon ii, consider the charged pseudo-scalar susceptibility \p S — J d 4 x(J^(x)J+(0)) (with J + = 
drf^u) as a function of /i/. This susceptibility can be represented as a functional integral; a straightforward computation 
using techniques similar to those used above, enables one to express this functional integral in terms of the spectrum 
of l0 (Ip+m): 

xUT,p I ) = yj ded^ (2e _^ (cogh(/3 (£ _ f ))+CQsm) =yj de (10) 

where the first equality is exact and the second equality exploits the fact that our interest is in the limit as T — > 0. 
From the second equality in eq. (|10|l we see that in the T — > limit, \p S diverges when 4r is increased to the smallest 
value of e for which the {p(e))o,ni is nonzero. From condition i) (previously demonstrated to hold) this implies that 
Xp S diverges when \xi equals the smallest value of e for which (/5(e))o,o is nonzero, namely e min . On the other hand, Xp S 
diverges in the infrared (at T = 0) when the lightest excitation with the quantum numbers of the ir + has an energy 
equal to fii. Clearly, if the state at T = and < were unchanged from the ordinary vacuum, this lightest 
excitation is simply m x . However, condition i) implies that in the absence of a first order transition, the system does 
remain in the vacuum state until {±j — 6 m i n . Together these imply that e m in = rn^ and condition ii) is established. 

In summary, the isospin Silver Blaze problem is resolved through the study of the spectrum of 70 times the usual 
Dirac operator. As T — > 0, the functional determinant for any gauge configuration is unchanged from its \i\ — value 
unless -^p exceeds the smallest value of |ej|, where €j is the real part of the eigenvalue of the operator. The resolution 
depends on the relevant gauge configurations at T — having a minimum \ greater than It was shown that 

this automatically happens in QCD unless there is a first-order phase transition. It would be interesting to see this 
picture emerging in lattice simulations via studies of the spectrum of 70 times the usual Dirac operator. 

Finally, as noted above there are analogous Silver Blaze problems associated with other chemical potentials. Of 
particular interest is the one associated with the baryon number: it is not plausible that we can understand nuclear 
matter in terms of QCD if we do not understand what happens below the phase transition to nuclear matter. The 
baryon Silver Blaze problem for small chemical potentials (i.e for hb < ^r^) is easily understood via the same methods 
as used here. However, the baryon Silver Blaze problem persists for all /is < Mn — B w 923MeV where Mm is the 
nucleon mass and B the binding energy per nucleon of nuclear matter. Moreover, in the regime, Mjv — B> [is > ^y*- , 
the functional determinants of the relevant configurations are altered from their \ij = values. Presumably, in this 
regime there are cancellations in the functional integrations due to the presence of the phase factors in eq. © which 
precisely compensate for the increases in the magnitude of the determinants. The understanding of how this comes 
about represents an essential challenge to resolving the baryon Silver Blaze problem and thereby gaining insight into 
the first-order phase transition to nuclear matter. 
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